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Abstract 

We calculate the finite temperature 2-loop effective potential in the MSSM with stop 
condensation, using a 3-dimensional effective theory. We find that in a part of the 
parameter space, a two-stage electroweak phase transition appears possible. The first 
stage would be the formation of a stop condensate, and the second stage is the transition 
to the standard electroweak minimum. The two-stage transition could significantly 
relax the baryon erasure bounds, but the parameter space allowing it {mn < 100 GeV, 
~ 155 — 160 GeV) is not very large. We estimate the reliability of our results 
using renormalization scale and gauge dependence. Finally we discuss some real-time 
aspects relevant for the viability of the two-stage scenario. 
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1 Introduction 



The electroweak (EW) phase transition in the early universe has been under intense 
investigation in recent years (for reviews, see [|l], 0). Indeed, the Standard Model (SM) 
offers the three necessary ingredients of Sakharov for baryogenesis However, 
careful studies have shown that in the SM (i) it may be somewhat difficult to produce 
enough baryons during the EW phase transition, due to the limited amount of CP- 
violation available and (ii) even if a 5 + L asymmetry is produced, it is washed out 
afterwards by the sphaleron effects. This is because the phase transition is only very 
weakly of first order for the experimentally allowed Higgs masses mn > 60 GeV, and 
even vanishes for mn ^rnw @]- 

Thus one has to return to speculative models of baryogenesis in the very early stages 
of the universe marked by grand unification and related topics — or one has to consider 
variants of the SM. The SM has been and still is incredibly successful in explaining 
experiments though its limited range is generally accepted. Supersymmetric exten- 
sions are today the only way to maintain perturbative calculability at higher scales 
and the practical success of the SM. Furthermore, from the theoretical point of view 
supersymmetry (SUSY) appears to be the only known attractive way of including quan- 
tum gravity through supersymmetric string theory. Still, SUSY has to be detected in 
experiments! 

The most simple SUSY model is the minimal supersymmetric extension (MSSM) 
of the SM: here the only new particles proposed are the superpartners with known 
(SM) couplings and an additional Higgs doublet. This is a very strong restriction in 
view of the multitude of possibilities offered, e.g., by string theory, but even in this 
simple case the necessary soft SUSY breaking terms introduce a huge number of free 
parameters. Fortunately we are interested in the behaviour of such a model at the 
high temperatures of a (modified) EW phase transition, in which case the problem 
simplifies. The new fermionic superpartners, the gauginos and higgsinos, become even 
more massive in the heat bath and can be integrated out. They may be important, 
however, for introducing strong CP-violating effects. The scalar quark partners, if 
heavy enough, can also be integrated out. Then the remaining nonzero Matsubara 
modes and heavy zero modes contribute to a 3-dimensional (3d) effective action like in 
the standard case [Q. "Integrating out" in the above means matching parameters in 
the comparison of a 4d and 3d set of amplitudes as in P, ^ §| . The effective theory for 
the light modes contains Higgs and gauge fields. If some combination of the two Higgs 
doublets is heavy compared with the scale g'^T, it can also be integrated out, resulting 
in an effective action of the same general form as for the SM [0, ^ |10|, This case 
leads to conclusions which are not much more optimistic than in the SM [9-14]: for 
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generic parameters, the phase transition would be at most relatively weakly of first 
order for Higgs masses rriH > 60 — 70 GeV. 

However, recently it has been argued that a rather light {mi^ < rrit) stop, the super- 
partner of the right-handed top, and low values of tan (3, allow for a strongly first-order 
phase transition even at Higgs masses rriH ~ rnw [0, |16|. Furthermore, the influ- 
ence of 2-loop contributions to the effective potential was observed to be large due 
to the strong interactions of the stop entering at the 2-loop level, strengthening the 



transition 17 



In this paper we investigate the MSSM with a light stop in some detail. In particular, 
we study the possibility of a two stage transition in which at some temperature the stop 
acquires a nonzero expectation value while the electroweak symmetry is still unbroken. 
In the second stage the universe would make a transition from the charge and colour 
breaking (CCB) minimum to the standard one. In Ref. [1^ it was argued that a 
CCB minimum does not develop during the electroweak phase transition as long as 
m^^ satisfies the bound necessary for vacuum stability at zero temperature; here we 
demonstrate that stop condensation may be possible in some part of the parameter 
space after all, leading to a two-stage transition in which the standard EW minimum 
is the global one at low temperatures. Such a transition would relax the baryon erasure 
bounds and allow for large Higgs masses, mn ^ 100 GeV. 

We compute the effective potential for the Higgs and stop fields at two loops in 
3d, working in an arbitrary background gauge. The gauge and renormalization scale 
dependence of our results should give an estimate of the size of higher loop contribu- 
tions. Finally, assuming that a two-stage transition takes place, we discuss some issues 
concerning the real-time history of such a transition. Previously, two-stage transitions 
(in two Higgs doublet models) have been considered in 



The plan of the paper is the following. In Sec. |] we discuss the 3d effective theory 
describing the EW phase transition in the MSSM. Sec. ^ contains the derivation of 
the 2-loop effective potential in the 3d theory. Implications from the 2-loop potential 
for the phase diagram of MSSM are presented in Sec. ^ In Sec. ^ we address some 
real-time aspects of the finite temperature phase transitions found, and the conclusions 
are in Sec. Some details related to Sees. ||, ^ are in the two appendices. 



2 The effective 3d theory and its parameters 

The effective 3d theory relevant for the EW phase transition in MSSM is |Tl[] 
L = \F:,F^, + \Gfpf^ + {DrH)\DrH) + m'^,H^H + \u,{H^Hf 

+ {DtU)\DtU)+ml^U^U + Xu:i{U^Uf + ^:iH^HU^U. (2.1) 
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Here = di - igw3r''A'^/2 and D| = di - are the SU(2) and SU(3) 

covariant derivatives, gw?, and 5^53 are the corresponding 3d gauge couphngs, H is 
the Higgs doublet and U is the right-handed stop field. This theory is an effective 
theory for the thermodynamics of the EW phase transition in MSSM, provided that the 
parameters are suitably fixed in terms of the temperature T and the zero-temperature 
physical parameters of the theory ||10| , pT| (see also below). The theory in should 
be particularly useful for the case that the right-handed stops are hght, mi^ < nit, 
which has recently attracted a lot of attention since the EW phase transition appears 



then to be strong enough for baryogenesis [T^, |T^, ||, |TT|. For a heavier stop, m^^ > rrit, 
the U field can be integrated out and the effective theory is just the 3d SU(2)+Higgs 
model as in the Standard Model P, 0. 

The reasons why it may be convenient (even in perturbation theory) to study the 
theory (|2.1| ) rather than to directly calculate the effective potential in 4d, are: 

1. The use of ( p.l|) factorizes the problem into two parts. The first part, the deriva- 
tion of (|2.1|) , is purely perturbative and free of IR-problems. The second part is the 
analysis of (|2.1|) in 3d and is subject to the usual IR-problems at finite T. Thus the 
IR-problems can be studied in a simplified setting. 

2. The effective theory (|2.1|) has a certain degree of universality, since the same 3d 
theory (with different parameter values) arises for many different 4d theories. Thus 
the IR-problems of many 4d theories can be studied once and for all using (|2.1| ). 

3. The construction of ( p.l|) automatically implements the daisy resummation pro- 
cedure required for the mass parameters at finite temperature, as well as resummations 
for all the couplings appearing in the theory (establishing the scale at which they are 
to be evaluated). In addition, there are resummations related to heavy modes in 3d. 
Apart from the daisy resummation, these resummations are usually not implemented 
in direct calculations in 4d. Yet their effect may be numerically non- negligible: e.g., 
the finite temperature strong gauge coupling is much smaller than gg{mz)T. 

4. One may eventually want to perform lattice simulations of the EW phase tran- 
sition in the MSSM by using (|2.1|) . Note that while the simulations in the Standard 
Model indicated that perturbation theory is surprisingly reliable for strong transi- 
tions [^, this may not necessarily hold in MSSM with two light scalar fields. 
In fact, we find that the renormalization scale and gauge dependence of the results is 
much larger than in the SM indicating that the convergence of the 2-loop perturbative 
expansion is worse. It should also be noted that there are always massless SU(3) gauge 
excitations related to an unbroken SU(2) subgroup in perturbation theory (even if the 
[/-field has an expectation value), unlike in the SM. If lattice simulations are ever made, 
then it is useful to have perturbative results in the same theory for comparison. 
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These benefits are slightly shadowed by the fact that at present there is no accurate 
derivation of the parameters in available for a wide range of 4d zero temperature 
parameters. While this derivation is straightforward and parallels the ones in P, [TO, 



11 1, it is complicated by the fact that the whole mass spectrum of the MSSM enters 



through relatively large radiative corrections, especially in the strongly interacting 
sector. Moreover (as is discussed in more detail in [Appendix A|) , the derivation should 
be made at 2-loop level as concerns the mass parameters, if one is studying a two-stage 
scenario where an accurate comparison of two critical temperatures is needed. This 
is because it is only a 2-loop calculation that fixes the scale appearing in the thermal 
screening terms (~ 93^^)- While a 2-loop derivation is quite doable 0, we will not go 
into it here. Instead we make a 1-loop derivation for a particular region of the parameter 
space, and an order of magnitude estimate of the 2-loop scale factors appearing. This 
derivation should illustrate with reasonable accuracy the real situation, and allows us 
to address some interesting problems: 



1. In |T^, it was found that a two-stage transition where the right-handed stop field 
{U) gets an expectation value should not be allowed, since then the CCB minimum 
would be the global one and one would remain there until the present day (a two-stage 
transition was argued to be ruled out in [^, as well, since it was found that there is no 
second order transition by which to get back to the standard EW minimum). Here we 
try to be somewhat more accurate in a particular region of the parameter space, and we 
demonstrate that a two-stage transition might be possible after all. While there remain 
uncertainties in our calculation and thus it is not clear whether this case is eventually 
realized, it is an interesting prospect and hence in our opinion worth studying. 



2. In [|T^ it was found that for small right-handed stop masses, the 2-loop corrections 
in the effective potential for the if-direction are large, making the transition stronger. 
Here we calculate the 2-loop potential in a general U, iJ-background. We confirm 
the 2-loop effect found in [|17|, find similar large 2-loop corrections in the t/-direction 
as well, and discuss the reliability of the statements based on these corrections. It 
should be stressed that the large 2-loop corrections are IR-sensitive 3d effects and thus 
independent of the 4d parameters used in the derivation of the effective theory. 

Since the derivation of the parameters in (|2.1| ) is a problem which completely fac- 
torizes from the calculation of the 2-loop effective potential in 3d, we will discuss this 
derivation separately in [Appendix A| . Let us here just summarize the zero temperature 
vacuum parameters used: the CP-odd Higgs mass is assumed for simplicity to be 
large (m^ ^300 GeV), the squark mixing parameters are very small, the running top 
mass is rrit ~ 170 GeV and the left-handed squark mass parameter is relatively small, 
niQ ~ 300 GeV. We then present the results as a function of tan/3 (or equivalently the 
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Higgs mass mn) and the right-handed stop mass m^-^. 

3 The 2-loop potential and physical observables 



The effective theory ( |2.1| ) contains two fundamental representation Higgs fields: the 
Higgs doublet H interacting via the electroweak SU(2) group, and the right-handed 
stop colour triplet U interacting via the strong SU(3) group. We have calculated the 
2-loop effective potential of this theory in a general H, U background. After gauge 
fixing, there remains a global symmetry in the theory, so the effective potential only 
depends on 

1 2 2 

H^H = ^, U^U=^. (3.1) 
2 ' 2 ^ ' 

In the 3d theory, the dimension of 0, x is GeV^/^ after a trivial rescaling with T; we 

nevertheless often use for clarity the 4d dimensions, so that 

04d _ 03d /„ r,N 

T " Vt' ^ ^ ^ 

and similarly for x- 

It should be noted that within the 3d theory ([2.1D, the temperature T does not 
appear explicitly. It is hidden in the parameters of the 3d theory, for instance in 
9s3 ~ ds'^- dimensionful 3d observables could hence be expressed in terms of 

(y^lg. However, to keep the connection to 4d physics clear and since the values of g^^, T 
are numerically close to each other, we will rather express the results in terms of T. 

We have calculated the effective potential ¥{(/), x) in a background field gauge with 
two gauge parameters, ^ related to SU(2) and ( related to SU(3). The applicability of 
the background field gauge to the present context has been discussed in . We display 



our results in the Landau gauge ^, C = 0, and other values of the gauge parameters are 
used in estimating the convergence of the calculation (see below). For generality, we 
have replaced SU(3) by SU(A^) in the calculation and we also keep a general dimension 
d in the formulas before the evaluation of the final two integrals ( [B.35| ), ( |B.40D , to which 



all others reduce. We work throughout in the MS scheme with the scale parameter /i. 
The details of the calculation are in [Appendix B. 



From the effective potential we extract several gauge fixing independent physical 
quantities. The critical temperature for the transition from the origin to the stan- 
dard electroweak (EW) minimum 0inin is defined such that at T^, 

V^(0,0) = V^(0^i„O). (3.3) 
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Similarly one finds the critical temperature for the transition from the origin to the 
charge and colour breaking (CCB) minimum Xmin in the U direction, and T^^'^ for the 
transition from the CCB minimum to the EW minimum. 

The strength of the phase transition with respect to the sphaleron rate is charac- 
terized by the vacuum expectation value (pmin/T in the broken phase. This is a gauge 
dependent quantity; we calculate it in the Landau gauge. A gauge independent char- 
acterization of the discontinuity could be obtained with the expectation values of the 
composite operators WH,WU [Q: 

d 



A{U^U) 



d 



^(0min,O)- 1^(0,0) 



dml^ L 



\/(0,0) 



(3.4) 



However, numerically A{WH) ~ 0Ln/2, A{WU) 
rately discuss the quantities in 
The latent heat is defined by 

d 



Xmin/2) SO that we do not sepa- 



L = T,— (p,(T)-pKT) 



T=Tc 



(3.5) 



where Ph{T) {pi{T)) is the pressure p{T) = — V^(min) in the high (low) temperature 
phase (note that we mostly write V{(j), x) in 4d units; the relation to 3d is V^d = TV^d, 
apart from field independent terms). For the surface tension we use the leading order 
expression 

a= / d^J2[V{v)-V{^,)], (3.6) 



where ip = {(f), x)y dip is the length element in the field space, and the path is the 
one that minimizes the result (this can be derived by extremizing the bounce action 
at Tc). Perturbative corrections to the surface tension from the derivative terms in the 
SU(2)-|-Higgs model were discussed in ||2^, p5| . 

It is important to estimate the reliability of the perturbative results obtained. To 
this end, one can use the fact that the physical observables derived are gauge- and 
/i-independent to the order to which they have been calculated. The remaining de- 
pendence then gives an estimate of the magnitude of higher-order corrections. For the 
gauge dependence, we vary ^, from zero to ~ 1. For the /2-dependence, we compare re- 
sults from the choice p. = T with results from a renormalization group (RG) improved 
optimized choice Q]. The optimization condition is that the 2-loop contribution to 
V'{ip) vanishes, and it gives a field dependent scale parameter fii^p) which should be 
proportional to some (non-linear) combination of the mass scales contributing to V{ip). 
The simple choice /i = T is also supposed to reflect the mass scales of the problem, 
since the non-perturbative 3d masses should be of order gf^g, g^^ [|, |2^, |2^, p8| . 
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Figure 1: The 2- loop effective potential at T = = = T^^'^ for tan/3 = 5 
{niH ~ 92 GeV), m^-^ = 158.3 GeV, /i = T, ^ = C = 0. 

4 Numerical results 

The 2-loop effective potential V{(1), x) for ^ particular choice of the parameters (tan/? = 
5, nii^ = 158.3 GeV, ft = T) at T = = 92.43 GeV is shown in Fig. This choice 
corresponds to a special case in which there is a simultaneous first order transition in 
both directions. In general, the transitions take place at different temperatures T/, T^, 
and there is a third temperature T^^'^ at which the broken minima are at equal heights. 

The 2-loop potential corresponding to the parameters of Fig. |1| is displayed more 
precisely in the 0- and x-directions in Fig. 0, together with the 1-loop potentials at the 
corresponding critical temperatures. One can observe several things. First, looking 
at the 1- and 2-loop potentials at the same /i = T, one sees that the 2-loop effects 
are large and seem to make the transition much stronger. For the ^-direction this 
was observed in |T^, and the reason for the strengthening was tracked down to 2-loop 



graphs involving the strongly interacting stops. From Fig. ^ one sees that there are 
similar large 2-loop corrections in the x-direction, as well. 

Second, comparing the 2-loop potentials at different choices of p, {p, = T vs. RG- 
improved /2), one can see that the reliability of the 2-loop results is questionable, 
especially for the 0-direction. The large /2-dependence arises from the /i-dependence 
of m^3(/2) (see Eq. ( |A.29| )) in the 1-loop corrections, and would be cancelled by 3- 



loop graphs involving strong interactions. Thus the 3-loop contributions are quite 
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Figure 2: The 1-loop and 2-loop effective potentials at the corresponding critical tem- 
peratures, for tan/3 = 5, m^^ = 158.3 GeV. It is seen that there may be a relatively large 
uncertainty in the prediction for the surface tension due to the large /i-dependence, es- 
pecially for the 0-direction. The figure has been drawn at the special point mi^ = mf^ ,, 
so that the critical temperatures T^~^ of the two transitions are equal (see Fig. |^). The 
/i-dependence is especially large at this point, since the mass parameters mjj^{fi), 
rriij^^fl) are close to zero so that the relative change can be significant. 
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important. The large /x-dependence seen here is in contrast to what has been observed 
in the pure SU(2)+Higgs model 0. At the same time, even in the SU(2)+Higgs 
model non-perturbative results for the surface tension (which is characterized by the 
size of the bump between the minima, see (|3.6| )) were significantly smaller than 2- loop 
perturbative results for weaker transitions 0, |2^ . Thus one can say that for the surface 
tension perturbative results are expected to give only an order of magnitude estimate. 
For the critical temperature the /2-dependence in Fig. |] is somewhat smaller, so the 
existence of the two transitions and their crossing point is on a firmer basis. For the 
vevs and latent heats the /i-dependence is also relatively smaller than for the surface 
tension, so the perturbative estimates should produce the correct qualitative features. 

The gauge dependence of V{(j), x) (obtained by varying ^, ( from zero to unity) is 
also shown in Fig. It differs from the /i-dependence in being much larger in the 
X-direction, especially for Tc. The reason for the difference is that the /i-dependence 
arises from scalar (in particular stop) excitations affecting also the 0-direction, whereas 
the gauge dependence arises from vector excitations. For quantities other than Tc the 
gauge dependence appears to be smaller. It should be noted that the gauge dependence 
(as well as the /i-dependence) is significantly smaller at the 2-loop than at the 1-loop 
level. The error bars containing an estimate of both the 2-loop gauge and /i-dependence 
for the critical temperatures are shown in Fig. ^. 

The dependence of the critical temperatures on m^^ and tan/5 is also shown in Fig. ^. 
The general pattern is the following: For a given tan/3, there is a certain critical value 
of m^^, for instance m^^ c = 158.3 GeV for tan/? = 5. If m^^ > rni^ ci the critical 
temperature of the normal (0 </>) EW phase transition is higher than the critical 
temperature of the CCB transition (0 — > x)- this case the EW phase transition 
proceeds in the normal way when the universe has cooled down to Tf. If m^^ < m^^ ^, 
on the contrary, then > and the universe ends up in the CCB minimum (provided 
that the first order phase transition is weak enough compared with the expansion rate 
of the universe). However, for the parameter values we are using, we find that the CCB 
minimum is not the global one at much lower temperatures. Instead, there is another 
first order phase transition at T^~*^ < in which the normal EW minimum becomes 
the global one. Thus, if this latter transition is weak enough to take place within the 
cosmological time scales available, one ends up finally in the normal EW minimum. 
The vacuum expectation value (j)rain can easily be large (pram/Tc ^ 1 after the latter 
transition (Fig. so that the sphaleron rate would be very effectively switched off. 

In Fig. ^ we have also shown the metastable branches of the EW (0 0) and CCB 
(0 x) transitions. These could be relevant if the transition with the higher Tc is so 
strong that it has not taken place before the temperature has cooled down to the Tc 
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Figure 3: The critical temperatures of the three transitions as a function of m^^ for 
tan/3 = 3, 12 (thin hues) and tan/3 = 5 (thick hues). The two-stage transition would 
take place to the left of the crossing point of the three critical curves; it is seen that 
a two-stage transition is possible, but there is not very much parameter space for it. 
The continuations of Tf to the left of the crossing point, and to the right of it, 
roughly represent the region of metastability of the transition with a higher Tc. Error 
bars indicate the dependences on fi,^,C (for tan/5 = 5), as explained in the text; the 
actual curves correspond to the Landau gauge and /2 = T. In terms of the parameter 
mfj = mf^ — m^, the x-axis is from —(70 GeV)^ to —(40 GeV)^. 

of the lower transition (see Sec. When one goes far enough into the metastability 
region, then the higher transition definitely does take place since it would have reached 
the barrier temperature by the time of the lower T^. 

5 Real-time history 

In the previous Section, we demonstrated that the phase diagram of the EW matter 
described by the MSSM may be such that the cosmological EW phase transition could 
take place in two stages. However, the transitions occurring would be of first order, so 
one has to discuss the amount of supercooling taking place and the possible reheating 
due to the latent heat released, as well, to see how the transitions would really proceed 
in the early universe. We will here make some simple estimates on the presumable 
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Figure 4: The expectation values 0min/^c, Xmin/^c of the three transitions (see Fig. ^ 
in the Landau gauge as a function of m^^ for tan/? = 3, 5, 12. The broken phase 
transition from Xmin > to 0min > is seen to easily lead to very large values of 0min- 
For clarity the metastable branches are not shown here. 

real-time events. 

The parameters needed for the real-time estimates are the latent heat L and the 
surface tension a, scaled with powers of Tc: L = L/T^, a = cr/Tj^ . Indeed, the small 



supercooling formula ^ for the nucleation temperature T„ = Tn/T^ is 

1-T„~0.4^, (5.1) 

assuming that in the high temperature phase p{T) oc T^. This formula can be expected 
to be valid only when 1 — T!„ <^ 1, and it usually breaks down by underestimating the 
value of 1 — T„. Another quantity of interest is the amount of heating that takes 
place after nucleation as the latent heat is released. Since the bubbles tend to fill the 
universe in a time scale very small compared with the time scale of expansion ^ , 
the reheating temperature Tj. can be estimated from 

e^iTn) = ei{Tr), (5.2) 

where eh(T),ei(T) are the energy densities in the high and low temperature phases, 
respectively. From (p.2|) one finds that reheating to the critical temperature (T^ > T^) 
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Figure 5: The latent heats of the transitions as a function of m^^ for tan/5 = 3, 5, 12. 
Note that at nii^^^, L^^'^ = L"^ — L^. For smaller rrif:^, L^^'^ grows rapidly. Metastable 
branches are not shown. 

takes place roughly ii L> 8cr^^^, assuming that the effective number of massless degrees 
of freedom is ~ 110. 

Results for the latent heat from ( ^.51 ) as a function of tan/?, m^^ are shown in Fig. ^. 

To compute the surface tension ( |3.6| ) one has to find the path in (</>, x) space which 
minimizes the integral. Near m^^ = fni^^^ the mass parameters of the 3d theory are 
small. Then the effective potential is dominated by the term 730^x^/4 in the tree level 
potential for 0, x 7^ (see Fig. |I|). Therefore the path which minimizes ( p. 61 ) for the 
transition x ^ </) is close to the rectangular path which first goes from (0, Xmin) to 
(0, 0) and then from (0, 0) to (0min, 0). We have used this path for estimating a''^^'^ in 
Fig. ^. The corresponding values are therefore an upper limit for the surface tension. 
At lower critical temperatures T^^"^, corresponding to smaller values of fn^^, the 3d 
mass parameters decrease and our simple approximation becomes less reliable. We 
have used different trial paths and compared the results for a^^*^ with those from the 
rectangular path. From this we estimate the error as at most 10% for the parameter 
range shown in Fig. ^. This error estimate, of course, does not contain the effects of 
higher-loop and non-perturbative contributions, which might make a^^'^ much smaller. 

Based on Figs. ^, ^ and on Eqs. ( |5.1| ), (|5.2|), one can say the following: 
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Figure 6: The surface tensions of the transitions as a function of m^^ for tan/3 = 
3, 5, 12. The curves for the broken phase x transition represent an upper bound 
obtained as explained in the text. Note that at m^^ ^, a^^'^ = cx'^ + a^. It is seen 
that for the broken phase transition, the surface tension becomes rapidly very large. 
Met ast able branches are not shown. 

1. The EW (0 0) and CCB (0 x) transitions seem to proceed in the normal 
manner: there is at most a few percent supercooling, and the latent heat released 
after the nucleation period reheats the matter considerably. Especially for larger Higgs 
masses in the EW phase transition, one might even reheat back to T/. 

2. For the broken phase transition {x ~^ the supercooling is considerably larger. 
For instance, at the point m^-^ ^ the latent heat is smaller [L^'^'^ = L'^ — L^) and the 
surface tension is larger {a^^'^ = a'^ + a^) than for the EW transition characterized by 
L'^, a'*'. Both effects increase the supercooling. The small supercooling formula ( [5.1| ) 
gives numbers of the order of 10%, but the true supercooling might be even larger. 
Moreover, when one goes to smaller the surface tension seems to grow more rapidly 
than the latent heat, so that one may end up in a situation where the lifetime of the 
metastable CCB minimum is comparable with cosmological time scales. This is clearly 



unacceptable. There is also no barrier temperature for the broken phase transition ^1 
Of course it might happen that some other nucleation mechanism enters, in which case 
the transition could proceed after all; or the non-perturbative surface tension might be 



13 



much smaller than our simple estimate. 

3. If one happens to be very close to m^^ ^, then things may become complicated. 
For instance, the supercooling for the CCB transition is generally somewhat larger 
than for the EW transition, so it is not excluded that at m^^ < m^^^^ one supercools 
so much that the metastable branch of in Fig. |^ becomes relevant. Then the EW 
transition might take place first. But after it has taken place, the CCB minimum is 
still the global one until T = T^^'f', so that in principle two other transitions (first to 
the CCB minimum and then back to the EW minimum at T < T^^^) could occur. 
However, this possibility clearly requires fine tuning. 



6 Discussion and conclusions 

We have studied the electroweak phase transition in MSSM with < rrit at the 
2-loop level of perturbation theory. It was found that in principle there is the pos- 
sibility that the EW phase transition takes place in two stages. The first part, the 
colour breaking transition, would delay and strengthen the transition to the standard 
electroweak minimum. The sphaleron rate should not be significantly switched off in 
the colour breaking minimum since the [/-field is an SU(2) singlet, but it would be very 
effectively switched off after the transition to the standard electroweak minimum for 
Higgs masses up to 100 GeV (in the standard minimum normal electroweak sphaleron 
estimates apply [0). Our calculations were performed in a particular region of the 
parameter space, with a relatively large CP-odd Higgs mass rriA and a relatively small 
left-handed third generation squark mass rriQ ~ 300 GeV, but it appears that relaxing 
these assumptions does not immediately change the qualitative pattern observed. 

The problem with the two stage scenario is that the parameter region in m^^ allowing 
for a two stage transition is rather small, m^^ ~ 155 — 160 GeV. There is always 
the danger that even though the CCB minimum becomes metastable, the first order 
transition to the standard EW minimum becomes so strong that the lifetime of the 
metastable state is comparable with cosmological time scales. Thus one may also 
interpret our result as a lower bound m^-^ ^155 GeV for the right-handed stop mass. 
This particular number is valid for small squark mixing parameters. 

Finally, it should be remembered that our study was based on 2-loop perturbation 
theory. In fact, already the large gauge and /i-dependence of our results indicates 
that higher-loop perturbative (and possibly non-perturbative) effects could be much 
larger than in the SU(2)-|-Higgs model which is relevant for the Standard Model and 
has been studied on the lattice. If the right-handed stop mass turns out to be such 
that the scenario studied could be realized, then one should probably consider lattice 
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simulations in the 3d SU(3)xSU(2) effective theory. 
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Appendix A 

In this Appendix, we discuss in some detail how the parameters in ( |2.1| ) can be fixed 
at 1-loop (and partly 2-loop) level. 



T = parameters 

Since our main point is studying 3d IR phenomena and since an accurate derivation 
of the 3d theory with the whole MSSM mass spectrum included would be somewhat 
complicated, we will consider a rather simple special case here. The final effective 
theory will be the same in many other cases as well, but the expressions for the 
parameters will be different. We consider the case of vanishing mixing in the squark 
mass matrix and a heavy CP-odd Higgs particle (m^ > 300 GeV) so that only one Higgs 
doublet remains; this is the most favourable case for baryogenesis . Neither of these 
assumptions affects the form of the final effective theory. For instance, for two Higgs 
doublets one can diagonalize the Higgs mass matrix close to Tc and integrate the heavy 
doublet out 0, ^, 0, |TT|. We mostly neglect the U(l) gauge coupling, which should 
be a good approximation with respect to the other uncertainties in the calculation. 
Finally, we assume that the only sparticles light enough to enter the calculation are 
the squarks in the third generation. Adding light fermionic degrees of freedom would 
again not change the form of the final effective theory, only its parameters (see [|TU| for 
the typical corrections arising). 

We shall for definiteness consider a scenario with a relatively light left-handed squark 
mass parameter mq ~ 300 GeV. This is at the lower end of the phenomenological 
constraints obtained with a realistic top mass JTB|, P^, and at the same time still allows 



the high-temperature expansion for T ~ 100 GeV (we also take the right-handed 
sbottom mass parameter to be ~ 300 GeV, although this has little effect). A 
much heavier mg does not change the pattern we are investigating qualitatively, but 
the formulas to be used simplify as the mg-corrections decouple. The mass parameter 
mfj is supposed to be relatively small, and its effects are discussed in the main text. 
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The most important renormalization effects of the T = vacuum theory with respect 
to our calculation concern the MS mass parameters m|^(/2), mfj{fi) (we work throughout 
in the MS scheme with the scale parameter /i). Calculating the pole Higgs mass rriH 
to 1-loop order, one finds that m'jj{p,) can be expressed in terms of physical quantities 
as 

mlifi) = + ^/.? sin^/5mj {^^^^ + ')^ (^.1) 



where only the dominant terms are kept. For instance, there is another term with 
rriQ — i> mfj and m^-^ ^tn^ ^^^^ much less important for the parameter values 
we are considering. For the relation between tan/3 and mn, we use the standard 1-loop 
formula containing the leading correction, 

m, 



4 = ml cos22/5 + In (A.2) 



Here, in the absence of mixing. 



2 2,2 2 2 , 2 , ^ 2 na /Ao\ 

^in =mu + mt, m^^ =mQ + mt+ -m^ cos2/5 . (A. 3) 

We take rrit = 170 GeV so that mf"^^ ~ 175 GeV according to m}'°'° ~ mt{l + ^f)- 
The bottom Yukawa coupling is neglected. 

Similarly, the dominant term in the expression for m^(/i) is 

mUfi) = ml + ^/i>j(21n-^ + l), (A.4) 

where mfj is expressed in terms of physical masses through ([A.3|) . 

Step 1: dimensional reduction 

Next we go to finite temperature. After integration over non-zero Matsubara frequen- 
cies, the effective theory for the thermodynamics of the electroweak phase transition 
contains only bosonic n = Matsubara modes [5-11,32]. The form of the effective 
theory can be written down immediately using 3d gauge invariance (only the terms 
coupled to H, U, A", Cf are shown explicitly below): 

L = \F^^F^^ + \Gfpf^ + {DrH)\DrH) + m],,H^H + \H,{H^Hf 

+ {DIU*)\DIU*) + ml^U^U + Xm{U^Uf + y^H^HU^U 

+ {DrQ)\DrQ) + {DtD*)\DtD*) + ml,QiQ^ + ml.DlD^ 

+ ^gwsQLQja sin^/3 HlHi + cos^P H^Hi {25iiSjk - SijSki) 
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1 2 



+ '^dsz^oUp DjD^ — Q*^Qjs (3(5q,5(5/3^ — Sa/sSjs) 

+ \gl,,H^HAlAl + \gl,C^Ci{U*)^\^\^U\ (A.5) 

where DfA^ = diA^ + gw'i^°'^^A\AQ and correspondingly for D^Cq, and H = iT2H*. 
We recall that, after trivial rescaling with T, the dimension of bosonic fields in 3d is 
QgYi/2 g^j^^ ^]-^g^^ Qf ^i^Q quartic couplings is GeV. In the interactions involving the fields 
Q,D,Aq, Cq which are to be integrated out shortly, we take the parameters only at 
tree-level. What then remains to be discussed are the values of the gauge couplings, 
the scalar couplings Xhs, A{/3,73, and the mass parameters. 

The gauge coupling in the dimensionally reduced SU(A^) gauge theory with Ns 
scalars and Nf fermions is [0, Q 



(A,6) 



where 



2. . _ 487r2 
^ ~ (22iV-iV5-4iVj)ln(/i/A)' ^ ^ 



and the standard logarithmic corrections from bosonic and fermionic integrals are 



Lkifi) = 21n^^21n- ^ 



47rT 7.0555T 



— 'y^ — 

Lf(a) = 21n^^21n ^ . (A.8) 

ttT 1.7639T ^ ' 

With the mass spectrum we are considering, for SU(3) the constants are = 3, 
Ns = Ng = 4 (four scalar triplets: U, D and the two components of Q), Nj = 6. For 
SU(2) we have N = 2, Ns = Ng = 4 (one Higgs doublet and the three components of 
Q), Nf = 6 (three families, each with four chiral doublets). 

The coupling g"^ is independent of /i to the order it has been calculated, but it 
depends on the temperature. To make this dependence simple, we rewrite ( |A.6| ) as 



9l = T 



'g', + ^^{22N-Ns-ANf)\n^ 



(A.9) 



where Tq = 100 GeV and gQ is g'^/T at the temperature Tq. For SU(3) with as{mz) ~ 
0.12 we get g^^ ~ 1.1, and for SU(2) with gwi^z) ~ 2/3 we get g^Q ~ 0.42. The effects 
of the U(l) gauge coupling are small and will hence be approximated by g' = 1/3. We 
denote §1^ = glJT, g"^^ = glrjT. 
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For the scalar self-coupling we only include the dominant /i^-correction, arising from 
the incomplete cancellation of bosonic and fermionic integrals at finite temperature: 



A 



H3 



= -(^7^ + (7'')cos22/3 +— /i^sinV ln2. (A.IO) 

For 73, A;73, it is more difficult to give the dominant radiative corrections. If one calcu- 
lates the first corrections with the particle spectrum we are using, then, for instance, 
one finds logarithmic terms in which are not the same as for the gauge coupling (7I3 
related directly to the gauge fields, although at tree- level \u3 oc g'g^. This is because 
supersymmetry has been broken by the particle masses and would be restored only at 
a higher scale. However, the relative effect of these corrections is smaller in 73, A[/3 
than in \h3i since the leading term is itself proportional to a large coupling. Hence it 
is sufficient for the present purpose to use the tree-level values 

Af/3 = ^^?|3, 73 = /i?sin'/3T, (A.ll) 

where /i^sin^/3 = g^m^ / {2m^) . 

The Higgs sector mass parameter m^j^ gets modified by the standard thermal screen- 
ing terms (e.g. fl^); and by a logarithmic Lb(/i)-term cancelling the running in ( [A.l|) : 



"i + ^/.?sin^/3m^f2 In 1:^ + 1 



2 167r^ ' ' 



L 



+ 



l-^{gl. + g")co,'2p +^^gl + l-^g'^ + \hlsin^p)T\ (A.12) 



The logarithmic running is numerically very significant, since it is proportional to the 
large parameter mq 3> m^. The scale at which the couplings in the screening part are 
evaluated can only be fixed with a 2-loop calculation, but since it is established by the 
same Lb(/2) and Lf{fi) functions as for the gauge coupling in ( [A. 6] ), we shall use the 
numerical values of g^^, g^r^ for g'siQw everywhere in the finite temperature formulas. 
Similarly, m\j^ is 

<3 = ml + ^^hlm%{2\n-^ + 1) + {^gl + ^hl{l + sin^P)^^ (A.13) 

Again the logarithmic running is numerically significant. 

For the other mass parameters we only include the leading screening terms: 

= ml+{^^gl + ^-gl + hi'l{l + sin^^)y\ (A.14) 
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m 



D3 



m: 



Co 



m 



2 

9w 



9l 



Nf m 



- H + 

3 6 



6 



Nf N: 



1 + — + — T' 



6 



6 



(A.15) 
(A.16) 
(A.17) 



where Nf = 6 (see the explanation below ([A.8|) ). Note, in particular, that there are 
similar logarithmic runnings in ( [A.14| ), ( [A.15| ) as in ( [A.12| ). However, for rriQ and mj^ 
these are not that important, as the parameters are themselves large. 



When is step 1 accurate? 

There are two basic requirements for the construction of the effective theory in ( [A.5|) 
to be accurate. First, the perturbative expansion for the parameters of the effective 
theory should converge. This expansion is free of IR-problems and thus proceeds just 
as perturbation theory at zero temperature, except that the mass scale is ~ 2ttT. Thus 
there should be no problems in the tan/3 region considered. Second, the higher-order 
operators neglected in (|A.5| ) should be insignificant. Such operators arise from the 
mass hierarchy of the scales kept and integrated out, in other words from the high- 
temperature expansion. Thus the effects are small provided that m^/(27rT)^ ^ 1, 
where m symbolizes the mass scales in the Lagrangian ( |A.5| ). This requirement is 
satisfied when the transition is not exceedingly strong. 



Step 2: heavy scale integrations 

To simplify ( |A.5| ), we will integrate out the 3d scales which are heavy at the transition 



rUl O. Afterwards, we make the replacement U* — U. 



point, namely Q, D, Ao,Co 
The result is the action in ( |2.1| ). 

The 1-loop relations between the parameters in and ( [A.5| ) are straightforward 
to derive using the techniques in [|^ ^ |T0|, |11[]. For the gauge couplings, one gets 



2 (new) 
9W3 



2 (new) 
^53 



2 

9W3 



9I3 



1 - 



'^9W3 

487rm^„ 

^3 



'^9W3 



487rmQ3 
Js3 



9I3 



487rm£i3 



(A.18) 
(A.19) 



Here the coefficient in the numerator of the first correction term is A^, and in the 
subsequent terms it is the number of components of the corresponding fields. 
The mass parameters and couplings change to be 



m 



2(new) 
H3 



m 



H3 



167r 



9w3''^Ao 



(A.20) 
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m 



2{new) 
U3 



A 



(new) 
H3 



^U3 - ^ah^Co - ^h'l3^Q3, 
^H3 - -T^T. 

3 1 



16 87rmQ3 



(7^3 cos^2/3 + Ag^^h^^ cos2/3 + 8hf^ sin' 



A 



(new) 
U3 



X 13_^|3 ^ 9S3 
■^U3 — 

36 87imcg 12 87rm£)3 



^53 



(new) 
73 



73 



h t3 



(A.21) 
, (A.22) 

(A.23) 
(A.24) 



When is step 2 accurate? 



The requirements for the integration in step 2 to be accurate are the same as for 
the dimensional reduction step. First, the perturbative expansion for the parameters 
should converge. The expansion parameters were estimated in 0], and are roughly 



9W3 



9I3 



for the integration over Aq, Cq, and 



^53 



7rmQ3 

for the integration over Q, D (for D, 



(A.25) 



(A.26) 

Evaluating the numerical values, 



one can see that when the 1-loop corrections are kept, the neglected 2-loop terms are 
small. Still, if the critical temperature is needed very precisely, the 2-loop corrections 
to the mass parameters would be needed (see below). 

The second requirement is that neglected higher-order operators should give small 
contributions. Such operators arise in particular from the expansion of the masses of 
the fields in ( |A.5| ) in the background (H'^H) = 0^/2, {U^U) = in terms of the 

background fields. Typical 1-loop terms to be expanded are 



1 



m: 



J3 + (4/^2 sin2/3 + gl, COS2/5 )0V8 + (S/i^ - g^x'/Q 
3/2 



3/2 



(A.27) 



The higher-order operators arise from the fourth terms in the expansions. In ( [A.27| ), 
the couplings and fields are for simplicity in 4d units. Taking the actual parameter 
values one can see that the expansion is valid for 0/T<3, x/T <3. 



20 



2-loop mass parameters 

So far we have discussed the mass parameters in ( |2.1| ) at the 1-loop level (~ g'gT'^). 
However, we are calculating the effective potential in 3d at the 2-loop level, and then 
the mass parameters get renormalized. In fact, the renormalized parts of the mass 
parameters in the MS scheme turn out to be of the form 

mlM = + Y^^Qa'ws + 9Ah3^7^3 - 12A^3 - 87! + 8^373) In (A.28) 

2 / — \ 2 / 4 \ 2 \ 2 2 2 \ 1 /a \ 

musW = + -[^[^933 + yAc73C/53 - 16A^3 - 273 + 3^1^373] In — -, (A.29) 
where m'jj,,m^, are the 1-loop expressions in (|A.2CI|) , ( |A.21| ). To fix Ah,, Af/3 would 



require a 2-loop derivation of the mass parameters (this is essentially a 2-loop calcu- 
lation of the effective potential in 4d; this is the only place where a 2-loop calculation 
in 4d gives information not contained in a 2-loop calculation in 3d [@, We will not 
make the 2-loop derivation here but shall instead take the order of magnitude estimate 

Ah3 ~ A{/3 ~ 7T. (A.30) 

This estimate arises from the typical mass scales of integrated-out degrees of freedom 
and can be confirmed in the case of the Standard Model, where a 2-loop derivation of 
the mass parameters has been made 0. The uncertainty from the choice of A's affects 
the critical temperatures, but does not affect the conclusions concerning the 2-loop 
effects on dimensionless ratios like 0/Tc, x/Tc, L = L/T^, a = u/T^-, where L is the 
latent heat and a is the surface tension. The estimated uncertainty in Tc is on the level 
of few (< 2 — 3) percent. 

This concludes the derivation of the parameters in the effective theory (p.l|). 



Appendix B 

In this Appendix, we present some details of the calculation of the 2-loop effective 
potential V{(j),x) in the theory defined by ( p.l| ) (with SU(3) replaced by SU(A^)). To 
compactify the expressions, we denote gwz — ^ Qw, 933 93, ^H3 ^h, ^U3 ^u- 

To derive V{(f),x), one has to shift the fields H and U by the constants H and U 
and then to calculate all the vacuum diagrams in the shifted theory |Q. Here we are 
using the notation 




For simplicity, we shall assume the shifts H and U to be real but otherwise arbitrary 
(the assumption of realness enters in the scalar propagators). We use the background 
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field gauge fixing conditions (the applicability of these has been discussed in p2|), 

G'^ = d^A'^ +'-^gw{H^r''H - HW'^h), (B.2) 
= d.Cf+'-(gs{u^\^U-U^\^u), (B.3) 

where r° are the Pauli matrices and A"^ are SU(A^) generators, normalized such that 
Tr A'^A'^ = 25"^^. The terms added to the action are 

^G^G^, (B.4) 

supplemented by the corresponding ghost terms. The limit ^ — 0, C 0, which we 
use for most of the numerical evaluations, gives the usual Landau gauge. 
The shift of H, U is made by writing the original fields as 

H = H + H, H, = ^{h. + ini), ^ = 1,2, 

U = U + JJ, u^ = a = l,...,N. (B.5) 

V2 

The gauge fields are not shifted. In the shifted theory, the mass spectrum is such that 
the (non-vanishing) gauge bosons masses are 

2 ^ 2 i2 2 -'-22 — 2 2(A' — 1) 2 /tj r\ 

= ■^9w(P . = l9sX , rn^ = — m^, (B.6) 

and the Goldstone boson masses are 

= + W + ^730' + C^G- (B.7) 

In the Higgs sector there is mixing, so that the eigenstates are 

1 



ml = - 
^ 2 



ml + ml± y'(m2 -m2)2 + 47|0V 



(B.8) 



where 



1 1 
1^1 = rriHs + ^^H<f + -73X^ ^1 = mfj^ + 3XuX^ + 2730^- (B.9) 



We denote the mixing angles by 



s^ = 4^^ c' = 4^^ cs = —?^. (B.IO) 
m+ — m_ m+ — m_ m+ — mi 
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The propagators in the shifted theory are as follows. Let 



Tap - Sap - -^^5 Lap - ^ (B-H) 



and correspondingly for the if-direction. Then 



- (B.13) 



For the vector and ghost propagators, we define the projection operators 

pAB ^ ^AB . l(N-2\ U^X^U U^X^U 1 U^Xi^Xm 
' 2\N-l) WU WU 2 f/tf/ ' ^ ■ ^ 

. + (B.18) 

pAB ^ ij^m^m^ 

2N-1 WU WU ' 
based on A^^^A:^^ = 2{6asSp'y - 6ap6^s/N) and satisfying ELi = ■ Then 

/ k k \ f P^^ P^^ P^^ 
{CfCf) = [5,j - — j {^-j^ ^ k^lml^ k^ + ml 

h h / r>^B pAB pAB s 



/ pAB pAB pAB X 

V k'' k^ + QrriQ k'^ + QmQj 

Here c^, are the ghost fields. For SU(2), P^^ = and m%,m% m^; Eqs. ( |B.20| ), 
( [B.21D then reduce to the standard isospin diagonal (5"'') propagators. 

The 1-loop potential of the theory can be easily calculated. Using P^"^ = 2{N — 1), 
Pg^^ = 1, the result is 



1 



m+ + m!_ + 3ml + 2(A^ - i)ml + ml 

B.22) 



127r 

+ 3(2 - e^')m'w + 2(iV - 1)(2 - e/')m% + (2 - e/')nrl 
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For the 2-loop potential, we need a number of integrals. Let 

S{e,m') = (B.23) 



Fijkip, q, r) = 6ik{pj - rj) + 5kj{ri - Qi) + 6ji{qk - Pk)- (B.25) 



Then we define the following integrals |35|, g, in terms of which the 2-loop potential 
can be written: 



Diy^iMlMlm') = 


1 dpdq {p\ Ml)V^^{q\ Ml)vUr\ Ml) 

y-Fijk{p, q, r)Firan{p, q, r), 
1 dpdq V;f(p2, M^)V}^{q\ Ml)S{r\ m^)A5,,5i^, 


(B.26) 
(B.27) 


Dgy{ml,mlM^) = 


j dpdqS{p\m\)S{q\ml)Vl^^{r\M^){p, - q,){p, - q,) 


,(B.28) 


D^QyimlmlM^) = 


jdpdqS{p\m\)S{q\ml)V^^{r\M^){-p,q,), 


(B.29) 


H{ml,ml,ml) = 


j dpdqS{p^,ml)S{q^,ml)S{r'^,ml), 


(B.30) 


4iM') = 


1 dpV^jip^ M')6,, = {d- l)Is{M') + (IsiCM^), 


(B.31) 


ls{m') = 


j dp S{p'^, m^), 


(B.32) 



where r = —p — g, is the space dimension, and 

All the D-integrals in ( |B.26| )-( |B.29| ) can be reduced to combinations of H- and Is- 
functions defined in ( |B.3CI| ), ( |B.32| ). In the Landau gauge, this reduction can be found 
in [^, and in the general gauge for mass combinations relevant for SU(2) in |2^ 
(note that in -Dyvv is —1/4 and -Dyvs is 1/4 of our corresponding expression). 
In the general gauge with mass combinations relevant for SU(A^), we have done the 
reduction using symbolic manipulation packages. We do not present the final formulas 
here since some of them are rather lengthy. As an example, let us give the most 
complicated Landau gauge result needed. For d = 3 — 2e, C = 0, 

+ H{,n^ M\ M^-) (-8M= - 10m= + 2^ + i^) 

+ H(m^,M^,0)i-2M^ + 5m^ -2— 

^ ' ' \ M2 2M4 2m2 
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m° A/f* 

+ ^(^^o,o)^ + i^(M^o,o)— 



where we have used that 



4M4 



I.(m^) = - — 



2m2 



m 



For the general gauge, let us give just the divergent parts of the /^-integrals: 



D'^y^iMlMlMl 



1 


1 




4e 


1 


1 


167r2 


4e 


1 


1 




4e 


1 


1 


167r2 


4e 



-(5 + 0{M',+M^ + Ml 
6 + 4C + 2C' , 
- 2m2 - 2ml 



These arise from the 3d if-integral [H defined in (p.30| ), 



Tj/ 2 2 2\ 



167r2 V4e 



h In 



mi + m2 + ma 2, 



(B.34) 
(B.35) 

(B.36) 
(B.37) 
(B.38) 
(B.39) 

(B.40) 



Using the Z)-integrals, one can write down the 2-loop graphs. Let us start with the 
nonabelian graphs, which we give for the SU(A^) sector; the SU(2) sector is a special 
case of this. There are the sunset graphs (PPP) and the figure-8 graphs (PP) (see, 
e.g., p2|, ^), containing vectors (P^V), scalars (P— *S) and ghosts (P^G): 

N 



(VVV) 



(GGV) 



-9l 



12 



(AT - l)(iV - 2)D^vv(0, 0, 0) + 3(iV - 2)D^vv('^G, ^g, 0) 



(B.41) 



(SSV) 



+3L'^vv("^G)"^G)"^g) 

[{N - l)(iV - 2)dSgv(0, 0, 0) + 2{N - 2)dSgv(0, C^^, ml) 
+ {N- 2)Di^^{Cm% (ml, 0) + 2D^^^y{Cml, (ml, ml) 
+Di^^{CmlXml,ml)^, (B.42) 

— D^ys(ml,ml,ml) 

^vvs(^G,^G,^')], (B.43) 

( ^2 _2n , /Ar iNnf 



(VVS) = -gl^ 



{N-l)Diy^{ml,ml,ml) + 2 



+N{N - 2)D^vs(0, ml, ml) + ^JL^ nC 



N 



-9l 



4 ^(iV - l)D^^^y{ml,ml,ml) + [N - l)/^LvK, ^g) 

~ -C 2 —2 \ I n'^ ^ — 2 _2 —2 



+NiN-2)Dgy{ml,ml,0) 



(BAA) 
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(GGS) 



(VV) 



(SV) 



2 2 2N-I 



H {(ml, (ml, ml) -H {(ml, (ml ml) 



N' 



m, 



%{CmlXml,ml) 



{N -1){N- 2)4(0)4(0) + 2{N - 2)/M0)/Mm^) 



l)I^{m'a)I^[maj 



9l 



G)'-vy"'G) + 

\\2N{N-2)li{^)Is{ml) 
+ {N- l)li{ml) [2Is{ml) + Is{ml) + /^(m^ 
+ [2/5(m^) + {N- l)Is{rnl) + (iV - l)Js(m^) 



(B.45) 



[B.46) 



[B.47) 



The integrals with all arguments zero vanish in the MS scheme (in a general scheme 
they give vacuum terms). We have also used the conventions 



/(..., m^,...) ^ s'f{. 



2 \ . 2 _[/ 2 
,m+, . . .) + c /(. . . ,m_ 

,<,...) + 



, m 



m_ 
2 



...)-/(•••,<,•••) 



(B.48) 
(B.49) 
(B.50) 



The last one is needed below. 

It is a useful check of (p.41 )-( B.47 ) to sum up the divergent parts using ( p.36|) - 
( P.39|) , and to verify that the result is gauge-independent. Indeed, the sum is 

AN^ -N + 3 



1 1 



{N' - 1) 



as^m + Y^573 + Y 



9l- 



4Ar2 



N +1 
N 



(B.51) 



The field independent part of this expression can be cancelled with a vacuum coun- 
terterm, and the rest is cancelled by the mass counterterms (see below). 

What remains is to write down the scalar graphs. The integrals appearing are sim- 
ple, but the result is complicated by the mixing of the two Higgs fields. To simplify the 
formulas, we employ ( [B.48| )-( |B.50| ). A scalar line involving the fields /i, vr is denoted 
by if, a line involving m, uj is denoted by f/, and a line involving the mixed propa- 
gator (p.l6| ) is denoted by X. We also leave out the subscript from the integral Is 
defined in (IB.32|). Then the results are: 



(UUU) 

(HHH) 
(HUU) 



-\lx^ [m{ml ml ml) + 2{N - l)H{ml ml, ml 

+H{ml,ml,ml), 
-A^02 [3H{ml, ml, ml) + 3H{ml, ml, ml) 

-^730' \Himl, ml, ml) + 2{N - l)H{ml, ml, ml) 



(B.52) 
(B.53) 
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+H (ml, ml, ml 



(UHH) 

(UXX) 

(HXX) 
(UUX) 

(HHX) 
(UHX) 

(XXX) 
(UU) 

(HH) 
(UH) 
(XX) 



- 2 2 

r73X 



H{ml, ml, ml) + 3H{ml, ml, ml] 



-{\^1<P'^ + ^^ul3X^)H{ml,ml,m 

^llx^ + 3^Hl3^^)H{ml,ml,ml), 

-Xul3<Px [3H{ml, ml, ml) + 2{N - l)H{ml, ml, ml) 

+H{ml,ml,ml) 
->^Hl3(l)x[3H{ml, ml, ml) + 3H{ml, ml, ml] 
-ll<PxH{ml,ml,ml), 



7I + 6XuXh) 4>xH{ml, ml, m: 



3Iiml)Iiml) + 4(iV - l)/(m2)/(m5) + 2/(m^)/(m2 ) 
+4iV(iV - l)I{ml)I{ml) + 4(iV - l)I{ml)Iiml) + 3I{ml)I{ml] 
3I{ml)I{ml) + 6I{ml)I{ml) + lhl{ml)l{ml] 
I{ml) + 2{N - l)I{ml) + I{ml)] [l{ml) + 31{ml) 
-f3l{ml)I{ml). 



1 

4 " 
1 



B.54) 
B.55) 
B.56) 
B.57) 

B.58) 
B.59) 
B.60) 

B.61) 

B.62) 
B.63) 
B.64) 
B.65) 



The complete 2-loop contribution to V{(j), x) is obtained by summing these scalar 
graphs together with the nonabelian graphs in (|B.41| )- (|B.47|) , evaluated both for SU(2) 
and SU(A^). 

The divergent part of the sum of the scalar graphs is 



(B.66) 



^ [A{N + l)Xtj + 273^ j + ^ ( UXj, + iV73 



When one sums this (for = 3) with (|B.51|) , evaluated both for SU(2) and SU(3), one 
gets exactly the divergences which after cancellation by the mass counterterms, leave 
for the renormalized mass parameters the /2-dependences seen in ([A.28| ), ( [A.29|) . 
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